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QUASILINEAR PARABOLIC STOCHASTIC PARTIAL DIFFERENTIAL 
EQUATIONS: EXISTENCE, UNIQUENESS 

MARTINA HOFMANOVA AND TUSHENG ZHANG 


Abstract. In this paper, we provide a direct approach to the existence and uniqueness 
of strong (in the probabilistic sense) and weak (in the PDE sense) solutions to quasilinear 
stochastic partial differential equations, which are neither monotone nor locally monotone. 


1. Introduction 

We consider a quasilinear parabolic stochastic partial differential equation of the form 
d u + di v(B(u)) d£ = di v(A(u)Vu) d t + a(u) d W(t), x G T d , t G [0, T], 

(L1) «(0)=uo, 

where W is a cylindrical Wiener process in H = L 2 (T d ) and a a mapping with values in the 
space of 7 -radonifying operators from H to certain Sobolev spaces. The coefficients B : R —> W 1 
and A : R —> M. dxd are nonlinear functions satisfying suitable hypotheses, in particular, the 
diffusion matrix A is uniformly elliptic. The precise description of the problem setting will be 
given in the next section. 

Equations of this type model the phenomenon of convection-diffusion of ideal fluids and 
therefore arise in a wide variety of important applications, including for instance two or three 
phase flows in porous media or sedimentation-consolidation processes (for a thorough exposition 
of this area given from a practical point of view we refer the reader to [GMT] and the references 
therein). The addition of a stochastic noise to this physical model is fully natural as it represents 
external perturbations or a lack of knowledge of certain physical parameters. 

Our aim is to establish existence of a unique solution to (1.1) that is strong in the probabilistic 
sense and weak in the PDEs sense. That is, we consider solutions that satisfy (1.1) with a given 
driving Wiener process and underlying stochastic basis in the sense of distributions. Recall that 
from the probabilistic point of view, two concepts of solution are typically considered in the 
theory of stochastic evolution equations, namely, pathwise (or strong) solutions and martingale 
(or weak) solutions. In the former notion the underlying probability space as well as the driving 
process is fixed in advance while in the latter case these stochastic elements become part of the 
solution of the problem. 

The existence and uniqueness of a weak pathwise solution, i.e. strong in the probabilistic 
sense and weak in the PDE sense, was obtained as a by-product in [DHV]. In this work, the 
authors were concerned with more general equations, namely degenerate parabolic SPDEs, and 
made use of the so-called kinetic approach and the notion of kinetic solution. The reason is 
that in such general cases, the degeneracy of the diffusion matrix introduces further difficulties 
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and the classical notion of PDE weak solution does not provide enough information in order to 
prove uniqueness and furthermore in some cases the equation might not even be well defined in 
the sense of distributions. The proof of existence in [DHV] relies on a Yamada-Watanabe-type 
argument (see e.g. [GK], [PR]). Pathwise uniqueness was established in the context of kinetic 
solutions which in particular implies uniqueness for PDE weak solutions in case they exist. 
Existence of a martingale solution was established in [DHV] via stochastic compactness method 
and these two results were combined and existence of a weak martingale solution deduced. 

In the present paper we put forward a direct (and therefore much simpler) approach towards 
existence and uniqueness of (1.1). First, we prove existence of a pathwise solution to (1.1). 
The proof is based on a suitable approximation procedure: existence of unique approximating 
solutions is established using the theory of locally monotone operators [LR] and then we show 
that these approximations satisfy several uniform bounds. In particular we prove that their 
gradients are uniformly bounded in space-time and have arbitrarily high moments. This part is 
based on the recent regularity result of [DDMH]. Finally, we are able to show strong convergence 
of the approximations and the limit is identified with a solution to (1.1). The proof of uniqueness 
is also new and much simpler than that in [DHV]. We believe the methods presented in this 
paper are also useful for tackling other type of fully nonlinear SPDEs. 

To conclude, let us mention several further references where similar problems were studied. 
In [DS] , analytical methods were used to prove existence and uniqueness of quasilinear parabolic 
SPDEs with second order operator having a dominating linear part and with gradient depen¬ 
dence in the noise. In the case of monotone coefficients, the literature is quite extensive, see for 
instance [PR], [BDPR1],[BDPR2] and the references therein. For locally monotone coefficients, 
see [LR] . For equations of gradient type see [Ge] . 

2. Mathematical framework 

2.1. Notations. In this paper, we adopt the following conventions. We work on a finite-time 
interval [0,T], T > 0, and consider periodic boundary conditions, that is, x € T d where T d = 
[0, l] d denotes the d-dimensional torus. Cl denotes the space of continuously differentiable 
functions, not necessarily bounded but having bounded first order derivative. For r £ [1, oo], L r 
are the Lebesgue spaces and the corresponding norm is denoted by || ■ ||ir. In order to measure 
higher regularity of functions (in the space variable) we make use of the Bessel potential spaces 
H a,r (T d ), net and r € (1, oo). Throughout the paper we will mostly work with the L 2 -scale 
and so we will write H a for H a ’ 2 (T d ) and H for L 2 (T d ). Recall that for all a > 0, the space H a 
is the usual Sobolev space of order a with the norm 



and that H a is the topological dual of H a . 

2.2. Hypotheses. Let us now introduce the precise setting of (1.1). We assume that the flux 
function 

B= (B 1 ,...,B d ) : M —> R d 

is of class Cl . The diffusion matrix A = {Aij) d J=l : R —>• R. dxd is of class C£, uniformly positive 
definite and bounded, i.e. <51 < A < Cl. 

Regarding the stochastic term, let (f2, , P) be a stochastic basis with a complete, right- 

continuous filtration. The driving process IT is a cylindrical Wiener process: it admits the 
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following decomposition 

OO 

( 2 . 2 ) W(t)=J2 

fe= 1 

where ( ek)k>i is some orthonormal system of the Hilbert space H , (/3k)k> l is a sequence of 
independent real-valued Brownian motions relative to (J^ t ). For each u £ H we consider a 
mapping a(u ) : H -7 H defined by a(u)ek = where CTfe(-) : R —> R are real-valued 

functions. In particular, we suppose that a satisfies the usual Lipschitz condition 

OO 

(2.3) \ a k(yi) - Vk(y 2 )\ 2 < C\yi - y 2 \ 2 . 

k= l 

This assumption implies in particular that cr maps H to L 2 (H, H) where L 2 {H 1 H) denotes the 
collection of Hilbert-Schmidt operators from H to H. Thus, given a predictable process u that 
belongs to L 2 (f2, L 2 (0, T; H)), the stochastic integral t i-7 J* a(u)dW is a well defined process 
taking values in H (see [DPZ, Chapter 4] for a thorough exposition). The equation (1.1) can be 
rewritten as 

OO 

dit + div(-B(w)) dt = div(kl(it)Vw) d t + <Jk{u) d f3k(t), x £ T d , t £ [0, T], 

( 2 - 4 ) ti 

u( 0 ) = uq. 


Later on it will be needed to ensure the existence of the stochastic integral in (1.1) as an 
7L a,r -valued process. We recall that the Bessel potential spaces H a ’ r with a > 0 and r £ [2, oo) 
belong to the class of 2-smooth Banach spaces and hence they are well suited for the stochastic 
Ito integration (see [Br], [BP] for the precise construction of the stochastic integral). So, let us 
denote by 7 (H, X) the space of the 7 -radonifying operators from H to a 2-smooth Banach space 
X. We recall that £ 7 (H,X) if the series 

^7fc^(efc) 

k> 0 


converges in L 2 (fl,X), for any sequence (jk)k>o of independent Gaussian real-valued random 
variables on a probability space (12, T, P) and any orthonormal basis {ek)k >0 of H. Then, this 
space is endowed with the norm 

'P £ j(H,X), 
xJ 


1*11 i{K,x) := E 


5^7 fc*(efc) 


k =1 


(which does not depend on ( 7 fe)fc>i, nor on (ek)k> 1 ) and is a Banach space. 

With this notation in hand, we state our last assumption upon the coefficient a. It coincides 
with the hypothesis (H a ^ r ) from [DDMH]: we assume that for all a < 2 and r £ [2, 00 ) 


(2.5) 


IK u )ll7 < 


C(l + IMIffo.r), 

C(l + ||u||ija,r + ||w||^l,ar), 


a £ [ 0 , 1 ], 
a > 1 . 


Detailed discussion of this condition was provided in [DDMH] so let us just make a few comments 
here. First of all we observe that there is an overlap between the Lipschitz assumption (2.3) 
and the assumption (2.5). For instance, it follows immediately that (2.3) implies (2.5) for a = 0 
and all r £ [2, 00 ). Nevertheless, for the purposes of our proof it proved useful to keep the two 
assumptions separate and thus we believe that it would not cause any confusion for the reader. 
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Indeed, assumption (2.3) is used several times throughout the paper whereas the use of (2.5) is 
somewhat hidden in Theorem 3.3 which is an application of the regularity result from [DDMH]. 
We end this section with a definition. 


Definition 2.1. An (J^ t )- adapted, H-valued continuous process ( u(t),t > 0) is said to be a 
solution to equation (2.4) if 

(i) u G L 2 (fl x [0, T], H 1 ) for any T > 0, 

(ii) for any (f> G C°°(Y d '), t > 0, the following holds almost surely 


( 2 . 6 ) 


{u{t),(j>) - {uoA) - f <B(u(s)),V$ds 
Jo 

= — f (z1(u(s))Vm(s), V<(>)ds + f (a(u(s)) dff(s), (j>). 
Jo Jo 


Remark that the solution is a weak solution in the sense of PDEs. 


3. Existence and uniqueness 

To begin with, write 

F(u) := — div(f?(w)) + div(A(M)Vu). 

We have the following estimate: 

(3.1) ll^Wllff-i <C(l + |M| ff i), ugh 1 . 

Indeed, for v G H 1 , it holds that 

|(F(u),n)| = | (B(u(x)),Vv(x)) - (A(ii(i))Vu(i), Vv(cc))| 

< c{ i + HfOiMiifi + c , ||u|| ff i||w|| ff i. 

This implies (3.1). 

Let P e ,e > 0 denote the semigroup on H generated by the Laplacian on T d . Recall that 

Pef{x)= [ P £ (x,z)f(z)dz , 

J T d 

here P E (x,z) stands for the heat kernel, x, z G T d . For p > 0, denote by C v (T d ) the space of 
functions that are ^-Holder continuous. We will use the following properties of the semigroup 
in the sequel. 

(3.2) \\P e f\\L-<C e \\f\\ H , f€H. 


(3.3) 



(Pe i (X, Z ) 


Pe 2 (x,z))h(z)dz 


< C'||^||c , ’)(T d )( e l 


S2) C 


h G C v (T d ), 


for some a v > 0. We refer the reader to [Gr] for these two properties. (3.3) can also be seen 
through the relation f Td P e (x, z)h(z) dz = E [h(Bf)\, where Bf is the Brownian motion on the 
torus T d . For e > 0, u G H, set 


A s (u)(x) = P e (A(u)){x), x G T d , 


(3.4) 

here 


P e {A(u)){x) = {P s {Ai j {u)){x))l j=1 . 
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Consider the following stochastic partial differential equation: 


(3.5) 


du e (f) + div(£?(u e (f))) d t = di v(A e (u e (t))Vu e (t)) dt + a(u e (t)) d W{t) 
u e { 0) = U 0 . 


Theorem 3.1. Let uq £ L 2 (fl,H). Then there exists a unique solution to the quasi-linear 
SPDE (3.5) that satisfies the following energy inequality: 

(3.6) sup /e sup ||u e (f)||^+ [ E||u e (f)||^i dtl < oo. 

e [ 0 <t<T Jo J 


Proof. First we claim that there exists a constant C such that 

(3.7) 5|^| 2 < A e {u){x)f, ■ £ < C|£| 2 for all e > 0, u £ H 1 , x £ T d , £ £ R d . 
By (H.2), one can find a constant C such that 

(3.8) <5|£| 2 < A(y)£-Z < C \£\ 2 for all y £ R, <£ e R d - 
Now, 

(3.9) A e (u)(x)£ ■ £ = [ P e (x,z)A(u(z))f ■ fdz. 

J T d 

Since f Td P e (x,z)dz = 1, (3.7) follows from (3.8) and (3.9). Set 

F E (u) := — div(B(u)) + div(T e (u)Vu), u £ H 1 . 

For u £ H 1 , by (3.7) we have 

{F E (u),u) = (7?(m), Vm) — (A e (u)'Vu,Vu) 

(3.10) <(C + C\\u\\ H \\u\\ m )-S\\u\\ 2 m 

<c+c\ 

for some constant <5i > 0. Moreover, 

(F s (u) -F e (v),u- v) 

= (B(u ) — B(v), V(u — v)) — (A e (m)Vm — ^4 e (u)Vw, V(m — v)) 

(3.11) = (B(u) - B(v), V (u - v)) - {A E (u)V(u - u), V(u - v)) 

- ((A e (u) — A e (v))Vv,S7(u — v)) 

<C + C\\u - v\\ 2 h - i5i||m - v\\ 2 h i - ((A e (u) - A E (v))X7v,V(u - v)). 
By (3.2) and the Lipschitz continuity of A we have 

~{(A e (u) - A e (v))Vv,V(u - v)) 

< \\Mu) - A s (m)||^oo( T d) ||m| j ||m - v\\ H i 

(3.12) < ||P £ [^(M)-T(M)]|| i cx ;(T d ) ||M|| H l||M-M|| ff l 

< Cell A(u) ~ A(v)\\ H \\v\\m\\ u ~ v \\m 

< C||m - v||/r||u||^i + (5 2 ||u — v\\%i, 

for some constant 62 < S±. 

Putting (3.11), (3.12) together we arrive at 

{F e (m) - F e (u),u-v) 

<C + C\\u — v\\~ H + C||m — u||^||m||^i — £> 3 ||u — u||^i, 


(3.13) 
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for some constant 63 > 0. (3.13) shows that F e satisfies the local monotonicity conditions 
imposed in [LR]. Applying Theorem 1.1 in [LR], we obtain the existence and uniqueness of the 
solution u £ . Next we prove the uniform bound in (3.6). By Ito’s formula, 

\\u e m 2 H = IKII 2 h - 2 f t (diy(B(u £ (s)),u £ (s))ds 
Jo 

+ 2 f (div(A £ (u £ (s))\/u £ (s)),u £ (s))ds 

Jo 

ft 00 

+ 2 (u £ (s),a(u e (s))dW(s)) + V / \\cr k (u £ {s) 

Jo Jo 


I H 


ds. 


fc =1 1 


By (3.10) it follows that 

\\u £ m 2 H +Si 


/Vooii 

Jo 


2 m ds 


(3.14) 


<IKII 2 H+ / (C + C||u e ( S )|| 2 ff )d S 

Jo 

rt 00 r t 


ft ft 

+ 2 / (u £ (s),a(u E (s))dW(s))+J2 ||cr fc (u e (s))||| f ds. 

Jo k=1 Jo 

Because the constants involved in the above equation are independent of e, the uniform bound 
(3.6) follows from (3.14), Burkholder’s inequality and Gronwall’s inequality. The proof is com¬ 
plete. □ 

Proposition 3.2. Let u 0 G LP(Q, T 0 ;LP(T d )) for some p G [2, 00 ). Then the solutions to (3.5) 
satisfy the following estimate 


(3.15) 


supE sup ||it e (t)||< 00 . 

£ 0 < t<T 


Proof. This result is obtained by a suitable version of the Ito formula using similar arguments 
as in the proof of (3.6). For further details we refer the reader to [DHV, Proposition 5.1]. □ 

As the next step we establish higher regularity of solutions to (3.5) which holds true uniformly 
in e. 


Theorem 3.3. Let uq G L p (Ll, Fq; C 1 +l (T d )) for some l > 0 and all p G [2, 00 ). Then it holds 
true that 

(3.16) supE sup ||Vu £ (t)||P oonrd , < 00 . 

Proof. The proof is based on Theorem 2.6 and Theorem 2.7 in [DDMH]. Since the setting in 
[DDMH] is slightly different, let us explain why the same ideas apply here. 

First of all, it is easy to observe that since the arguments for Dirichlet boundary conditions 
are more involved, considering periodic boundary conditions simplifies the proofs and does not 
cause any additional difficulties. The main difference between our equation (3.5) and the model 
problem from [DDMH] is that our second order operator A e is by definition nonlocal. Let us 
thus repeat the main ideas from [DDMH] and justify each step. 

We consider the following auxiliary problem 

dz £ = A z £ di + <r(u e ) dW, 

* e (o) = o, 
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and define y E = u e — z E . Than y e solves 

d t y e = div(A e (M e )Vy e ) + dw(B{u E )) + div((H e (tt e ) - I)Vz e ), 

2/(0) = Mo, 


which is a (pathwise) deterministic linear parabolic PDE. 

To establish the first step in the regularity problem, i.e. [DDMH, Theorem 2.6], we remark 
that all the corresponding estimates for the stochastic convolution z E are valid uniformly in £ due 
to (3.6) and (3.15). Moreover, the estimates for y E depend on A e {u E ) only through the ellipticity 
and boundedness constants from (3.7) and therefore is also independent of e. Consequently, we 
deduce that there exists r\ £ (0,1) such that for all p £ [2, oo) 

(3.17) supE||M £ ||^ ([0T]><Td) <oo. 


We proceed with the next step proven in [DDMH, Theorem 2.7, case k = 1]. The same 
arguments as above apply to the bounds of the stochastic convolution here. However, concerning 
the estimates of y E one has to be more careful. In view of [DDMH, Theorem 3.3] we need to 
verify that 

(3.18) A e {u E ), B{u E ), ( A e (u e ) - I)Vz e G 77(12; C a/2 ’ a ([ 0, T\ x T d )) 

for some a £ (0,1), where C“/ 2, “([0, T] x T d ) denotes the space of functions that are a-Holder 
continuous with respect to the parabolic distance 

d((t, x ), (s, y)) = max{|t - s| 1/2 , \x - y |}. 


To this end, we observe that if f £ C a//2, “([0, T] x T d ) then P e f £ C“/ 2,a ([0, T] x T d ) uniformly 
in £. Indeed, since the convolution kernel P e ( x, z ) depends only on the difference x — z, we have 


\P £ f(t,x) - P s f(s,y)\ _ | f Td P e (x - z)f(t,z)dz- f Jd Pe(y - z)f(s,z)dz\ 


d((t,x), (s,y)) 


< [ Pe{z) 
JT d 


d((t,x), ( s,y )) 

| f(t,x-z) - f{s,y — z)\ 
d((t,x), ( s,y )) 


dz 


7 ll/llc“/ 2 .“([0,T]xT< J )- 


Therefore due to (3.17) and the Lipschitz continuity of A we conclude that A e (u E ) as well as 
(. A e {u E ) — I)Vz £ possess the regularity required in (3.18) uniformly in e. The corresponding 
statement for B(u e ) follows immediately since B is Lipschitz. Finally, [DDMH, Theorem 2.7] 
applies in particular (3.16) follows. □ 


Theorem 3.4. Let uq £ L m (Q, To; C 1+ *(T d )) for some l > 0 and all m G [2,oo). Then 
there exists a unique solution to the quasi-linear SPDE (1.1) that satisfies the following energy 
inequality 

(3.19) E sup \\u(t)\\ 2 H + [ E||u(t)||^-i dt < oo. 

0 <t<T Jo 


Proof. We first establish the existence. Let u E be the solution to equation (3.5). We will show 
that u E converges to a solution to equation (2.4). The estimate (3.6) implies that there exist a 
sequence ( u En ,n >1), and a process 

u £ l 2 (vl x [o,ir], jt 1 ) n i 2 (f2,^^(o,r, jt)), 

for which the following holds: 

(i) u En —> u weakly in L 2 (fl x [0,T], H 1 ), hence weakly in L 2 (LL x [0,T],i7). 
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(ii) u En —» u in L 2 (fl, L°°(0, T, H)) with respect to the weak star topology, 

Next we show that u E actually converges to u in L 1 (n,i?) as e —>• 0. It is sufficient to prove 
that u En is a Cauchy sequence. 

Let 0 < £1 < £ 2 - By Ito’s formula, 


(3.20) 


r t 


= 2 f (B(u £ 1 (s)) — B(u £ 2 (s)),V(u £ 1 (s) — u £ 2 (s)))ds 

J o 

— 2 f (A £l (u £l (s))'Vu £l (s) — A £ 2 (u £ 2 (s))\7u £ 2 (s),\7(u £l (s) — u £ 2 (s)))ds 
Jo 

+ 2 f (u £l (s) — u £2 (s), (a(u £l (s)) — a(u £2 (s))) dW(s)) 

Jo 


f 0 

oo »t 

+ E / ||a fe (^(«))-^(« £2 (s))|||d S 

fc=l "'° 

:= I\(t) + /2(f) + 13(f) + /4(f). 


By the Lipschitz continuity of f?, for any <5i > 0 there exists a constant C\ such that 


(3.21) 


/i(f) <Ci [ ||u £l (s) - u £2 (s)llf^ds + <5i f 
Jo Jo 


|u £l (s) 


w £2 (s) lipids. 


By (3.7), we have 


(3.22) 


d 2 (f) = —2 [ (A El (u E1 (s))V(u ei (s) — u e2 (s)), V(u £1 (s) — u £2 (s)))ds 
J 0 

-2 [ ((A ei (u El (s)) - A £2 (u £2 (s)))V« £2 (s), V(u £l (s) - u £2 (s)))ds 

Jo 

<—2<5 f ||it £l (s) — tt £2 (s)lipids 

Jo 

- 2 [\(A E1 (u £l (a)) - A £2 (u £l ( S )))V« £2 (a), V(u £l (a) - u £2 (a))>da 

Jo 

- 2 [*((A E2 (u ei (a)) - 4(u £2 ( S )))V^O),V(« £ ‘( S ) - u £2 (s)))d S . 

Jo 
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Let S 2 be a small constant to be fixed later. In view of (3.3) we have 


(3.23) 


-2 [ (( A eAu ei (s)) - A Ea (u ei (s)))Vu ea (s), V(u ei (s) - u £2 (s)))ds 
Jo 

<S 2 [ ||w El (s)-w £2 (s)||^ids 
Jo 

+ C [ f \A Sl (u Sl (s))(x) — A e2 (u £l (s))(a;)| 2 1 Vu E2 (s)| 2 (x) dxds 
Jo Jf d 

<S 2 [ ||m Ei (s) - w E2 (s)||^ids 
Jo 


|Vzx E 2 (s)| 2 (x) drrds 


+ C f f f (P ei (x,z) - P e 2 (x,z))A(u ei {s))(z)dz 
Jo Jf d JT d 

<8 2 f ||n £l (s) — w £2 (s)||^ids 

Jo 

+ C J J d \\A{u ei {s))\\ 2 Cri( j d) {e\ - e 2 ) 2 “”|Vu £2 (s)| 2 (a;)da;ds 

<S 2 ( ||w £l (s) - M £2 (s)||^ids + C'(ei -e 2 f olr ' [ ||w £2 (s)||^i (l + ||u £l (s)||^ (Td) )ds. 
Jo Jo 


Due to Lipschitz continuity of A 1 for any positive constant S 3 we have 


(3.24) 


-2 [ ((A e 2 {u ei (s)) - A £2 (« £2 (s)))Vit £2 (s), V(u £l (s) - u £2 (s)))ds 
J 0 

<83 f ||u £l (s)-u £ 2 (s)||^ids 
Jo 

+ C [ [ \A E 2 (u ei (s)){x) - A £2 (u £2 (s))(x)| 2 |Vu £2 (s)| 2 (ai)da:ds 
Jo J T d 

<83 [ ||u £l (s)-it £2 (s)lipids 
Jo 

+ C f ||Vu £2 (s)|||oo f \P S2 [A(u^(s)) - A(n £2 ( s ))](x)| 2 dxd S 
Jo Jj d 

<83 [ ||u £l (s)-u £ 2 (s)||^ids 
Jo 

+ C j 1 ||Vu £ 2 (s)|| 2 tc(t < 1 ) ||u £i (s) - u £ 2 (s)||^ds. 

J 0 
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Choose <5 2 , 63 sufficiently small and substitute (3.24), (3.23) into (3.22) to obtain 

hit) < -64 f ||u £1 (s) - u £2 (s)||^ids 
Jo 

(3-25) + C( £l - e 2 ) 2 ^ f \\u E2 (s)\\ 2 H1 (l + ll^llc^))^ 

Jo 

+ C /* t ;||Vu £ 2 (s)|| 2 00 ||u £l (s) — u £ 2 (s)||| r ds, 

Jo 

for some (5 4 > 0. Choosing <5i < 64 it follows from (3.20), (3.21) and (3.25) that 

IIu £l (t) - < -(«4 - s 1 ) f ||u £l (s) - U £2 (s)||^ds 

Jo 

+ C [ ||u £l (s) -U £2 (s)||^ds 

Jo 


(3-26) +C( £l -e 2 ) 2 ^ f (1 + ||u £l WII^ (T<J ))h £2 ( S )||^d S 

Jo 

+ C [\\Vu E2 (s)\\ 2 L „\\u E '(s) -u E2 (s)\\ 2 H ds 
Jo 

+ 2 f (u £l (s) — u £2 (s), (ct(u £i (s)) — cr(u £2 (s)))dVF(s)). 

Jo 


For any M > 0, define 

t m 2 = inf {* > 0; ||Vu £2 (t)||i= > M or ||u £l (t)||c»>(T<*) > M} 

with the convention inf 0 = T. Then is an (J^-stopping time. Keeping the bound (3.6) 
mind and replacing t by t A t]j 2 in (3.26) we deduce that 

E||u £l (t A rlf) — u E2 {t A t\' 2 )\\ 2 h 


(3.27) 


< C [ E||u £l (s At^ 2 ) -u £2 (sAr^ 2 )||^ds 
Jo 

r tAr M 2 

+ C M E / ||u £2 (s)||^ids(ei - e 2 ) 2a,) 


+ CM 2 I E\\u ei (s Arif)-u e2 (s Arif )\\ 2 H ds. 


JO 

By the Gronwall’s inequality we obtain from (3.27) that 

(3.28) E||it £l (f A t^ 2 ) — u E2 (t A t]^)\\ 2 h < Cm{si — £ 2 ) 2a '' exp(CT + CM 2 T). 


By (3.6), 

E\\u E '(t)-u E2 (t)\\ H 

= E[||u £l (t A r If) - u E2 (t A T^ 2 )||l r i, 2 > t ] + E[||u £l (t) - u E2 (t)\\ H l T ^ <t ] 

< E||u £l (t A rlf) -u E2 {t At^ 2 )|| h + C'sup(E||u e (t)|||-)5[P(T]) 2 < t)]’ 

£ 

<E\\u El {t A rl; 2 ) - u E2 {t A t]j 2 )\\h + C[E{rlf <t)]i. 


(3.29) 
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Given any rj > 0. In view of (3.16), (3.17) and (3.6) we can first choose M such that 
C[P(t],} 2 <# 

<C{[P(sup ||Vu £ 2 (s)|| I ,«. (T i ) > M )]5 + P( sup ||M ei (s)||c'»(T‘i) > M)]?} 

0 <S<T 0 <S<T 

77 

< —, for all e 2 , £± > 0 . 

Then use (3.28) to find £o so that for £±,£2 < £o? 

(3.30) E\\u ei (t Arif)-u e2 {t/\rlf)\\ H < |. 

Because 77 is arbitrary, we conclude from (3.29), (3.30) and (3.30) that for all t £ [0, T], u e (t) —> 
u(t) in L 1 (0, H) and according to (3.6) and Vitali’s convergence theorem, we deduce that u E —» u 
in [0, T], H), which can be further improved using (3.15) and Vitali’s convergence theorem 

to 


(3.31) 


u E —} u in L P (Q x [0,T], H) VpE[l,oo). 


Next we show that the limit process it is a solution to equation (2.4). To this end, take a test 
function (j> £ C°°{ T d ) and use equation (3.5) to get 


(3.32) 


(u e (t), (j>) - (u 0 , <j>) - [ (B(u E {s )), V</>)ds 
Jo 

= - f (A e (u E (s))Vu E (s), V<£)ds + [ (a(u E (s))dW(s),(p). 
Jo Jo 


Taking (3.31) into account and letting e —> 0 in (3.32) we obtain 
E|(« e (t) — u(t), 4>)\ —> 0, 


E 


( B(u e ) — B(u),V(/))ds 


< IIV^HiocE f || u e — u||_f/ds —* 0 . 
Jo 


For the stochastic integral, we have 


E 


r>t f ft 0 o 

/ ((a(u e ) - a(u))dW,(/)) < CE( / ^(a k (u E ) - a k (u), (j>) 2 ds 
Jo \J 0 k=1 


<CU\\ H E^j \\u E -uf H ds 

It remains to pass to the limit in the second order term. Write 

f (A e (u e (s))Vw E (s) — A(u(s))Vu(s), V</>)ds 
Jo 

= [ {( A e{u e (s)) - A e (u{s)))S7u e {s), V</>)ds 
Jo 

+ [ ((Ae(u(s)) - A(u{s)))Vu E (s), V</>)ds 
Jo 

+ f (A(u(s))(Vu E (s) — Vu(s)), V<fi)ds. 

Jo 


0 . 


(3.33) 
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By the contraction property of the semigroup P E , Lipschitz continuity of A , (3.6) and (3.31), we 
have 


E 


((A E (u e (s)) - 4(«(s)))V« £ (s), S7</>)ds 


(3.34) < ||V</)||z,=oE [ f \P e (A(u e (s)) - ^(u(s)))(a;)||Vu e (s)(a;)| dzds 

Jo Jt* 


<C(EjJ\A(u%s))-A(u(s))r H ds) ^J o nu e (s)\\ 2 Hl dsj -M). 

By the strong continuity of the semigroup P e and boundedness of A(u), we have 


E 


(( A e (u(s )) - 4(it(s)))Vti e (s), V</>)ds 


(3.35) < ||V0|| L =o( T d)E / f \P e (A(u(s))(x) - A(u(s))(a;)||Vu £ (s)(j;)| dxds 

Jo JT d 


<C(E J \\P e A(u(a)) - A(u(s))r H ds 


,£ ( s )llff 1 ds ) 0 . 


By the weak convergence of u e , 

(3.36) E [ (A(u(s))(Vit e (s)-Vu(s)),V0)ds^O. 

Jo 

Putting together (3.33)-(3.36) we arrive at 

(u(t), <j>) - (u 0 , 4>) - f ( B(u(s )), V</>)ds 
Jo 

= — f (A(u(s))Vu(s), V</>)ds + f (a(u(s))dW(s),<t>) 

Jo Jo 

proving the existence of a solution. 

Next we prove the uniqueness. Let 1 > a\ > a ,2 > • • • > a n ■ ■ ■ > 0 be a fixed sequence of 
decreasing positive numbers such that 

r i 


— du = 1, • • • 


— d u = n, 
u 


Let i/j n (u) be a continuous function such that support C (a n ,a n - 1 ) and 

^ ^ rCln — 1 

0 < 4>n(u) < 2- X 

n u 


Define 

We have 
(3.37) 


/ ip n (u)du = 1 . 

a n 

/•Pi rv 

<t>n{x) = / / ipn{u)dudy. 

Jo Jo 


\<l>n( x )\ < !. 0 < 0"(x) < 2- X -i-. 

n x 
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Introduce a functional <h n : H —>• K by 


Then we have 


and 


(3.38) 


$ 


i(u) = / <p n (u(z)) dz, 

JT d 


u€ H. 


K( u )( h )= [ <t>n(u{z))h{z)dz, 

Jf d 


KMMd) = [ <f>n(.u(z))h(z)g(z)dz. 
JT d 


Suppose that ui,u 2 are two solutions to equation (2.4). We may apply the generalized Ito 
formula [DHV, Proposition A.l] to deduce 


(3.39) 


$n(ui(t) - u 2 (t)) = / </>„(wi(f, z) -u 2 (t, z))dz 

Jf d 

= / $(j(ui(s) -u 2 (s))(-div(S(ui(s))) + div(_B(u 2 (s))))ds 
J o 

+ / $^(ui(s) — u 2 (s))(div(A(ui(s))Vui(s)) — div(A(it 2 (s))Vw 2 (s)))ds 
Jo 

+ [ $n( u i( s ) - w 2 (s))(ct(ui(s)) -0-(tii(s)))dW(s) 

Jo 

+ \ J Q tr [( Cr ( U l( S )) - Cr ( M l( S )))* ° $ "( U l( S ) - Ms)) ° (ct(mi(s)) - CT(ui(s)))]ds 

= I 1 n (t) + In(t)+Mt)+Mt). 


We will bound each of the terms on the right. Now 


(3.40) 

In(t) = / / <t>n{M s , z ) - u 2 (s,z))V(ui(s, z) - u 2 {s,z)) ■ B(ui(s, z)) - B(u 2 (s, z)) dzds 

Jo JT d 

<— f [ |V(ui(s, z) -u 2 (s,z))|d.zds 
n Jo JT d 
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where the Lipschitz continuity of B and (3.37) have been used. For I 2 , we have 

I n(t) = ~[ [ <t>n(ui(s,z) -U 2 {s,z))V(ui(s,z) - U 2 (s, z)) 

Jo JT d 

■ A{u\(s, z))V(ui(s, z ) — U 2 {s, z)) dzds 


(3.41) 


K(ui(s,z) -u 2 (s,z))V(u 1 (s,z) - u 2 (s,z )) 

/ 0 J T d 

■ (A(ui(s, z)) - A(u 2 (s, z)))Vu 2 (s, z) dzds 

r t 


< -6 


[ [ 0n{u 1 (s,z)-u 2 {s,z))\v(u 1 (s,z)-u 2 (s,z))\ 2 dzds 
Jo j T d 

+ — [ [ |V(ui(s,2;) - u 2 (s, z))| |Vu 2 (s, ^)|dzds 
n Jo J T d 

<s[ [ <j>'^{ui{s, z) — u 2 {s, z))|V(mi(s, z) — u 2 (s, z))\ 2 dzds 
Jo j T d 

+ — [ [ |Vu 2 (s, z)\ 2 d~ds + — / / |V«i(s, z)\ 2 dzds, 

n Jo Jt 1 n Jo Jj d 

where the Lipschitz continuity of A and (3.37) have been used. The fourth term in (3.38) can 
be estimated as follows. 

oo 

^nO) = ( $ "(wi(s) - U 2 {s)) O (ct(mi(s)) - a(u!(s)))e k , (ct(ui(s)) - o-(« 1 (a)))e*)ds 

Jo 


k= 1 
oo r t 


=f-r r 


(3.42) 


k=1 Jo Jf d 


o Jf d 


</)"(iti(s, z) - u 2 (s, z))l(a(u 1 (s)) - tr(Mi(s)))e fc ( 2 )| 2 dzds 

OO 

</)”(ui(s,z) -m 2 (s,z))^ |(<7fc(ui(s)) - a k (ui(s)))e k (z)\ 2 dzds 


fc=l 


<— [ [ |mi(s, z) — u 2 (s, z)\ dzds. 
n Jo J T d 

Substituting (3.40), (3.41) and (3.42) into (3.39) we get 


E4>„(ui(f) - u 2 (t)) = E f </>„(ui(i, z) - u 2 (t, z)) dz 
Jf d 

‘ f [ |Vui(s, z)\ 2 dzds + E f ( |Vw 2 (s, z)\ 2 dzds 
Jo Jr 1 Jo J T d J 

) / |ui(s, z) — u 2 (s, z )\dzds 

Jo Jt 1 


< —I E 

n 


Hence by (3.19), 


Letting n —> oo we obtain 


r C 

E cj) n (ui(t, z) -u 2 {t, z))dz < —. 
Jfd n 


e[ \ui(t,z) — u 2 (t,z)\dz = 0 
Ji d 

This completes the proof of the theorem. 


□ 
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